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We propose a gravimeter based on a matter-wave resonant cavity loaded with a Bose-Einstein 
condensate and closed with a sequence of periodic Raman pulses. The gravimeter sensitivity in- 
creases quickly with the number of cycles experienced by the condensate inside the cavity. The 
matter wave is refocused thanks to a spherical wave-front of the Raman pulses. This provides a 
transverse confinement of the condensate which is discussed in terms of a stability analysis. We 
develop the analogy of this device with a resonator in momentum space for matter waves. 
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The realization of gravimeters loaded with cold atomic 
clouds has drastically increased the accuracy of the mea- 
surement of gravitational acceleration to a few parts per 
billion [U [2]. The recent obtention of quasi-continuous 
atom lasers now opens new perspectives for gravito- 
inertial sensors with the possibility to load these devices 
with a fully coherent and collimated matter source in- 
stead of the incoherent cold atomic samples used so far. 
In this paper we investigate a gravimeter based on a res- 
onant matter-wave cavity loaded with a Bose-Einstein 
condensate. The condensate is stabilized in momentum 
space thanks to a sequence of periodic "mirror pulses" 
consisting in velocity-sensitive double Raman 7r-pulses. 
Between the pulses, the optical potential is shut off and 
the condensate experiences a pure free fall. Other mea- 
surements of the acceleration of gravity have been pro- 
posed, based on the Bloch oscillations of an atomic cloud 
in a standing light wave [SI H] or on the bouncing of a 
cloud on an evanescent- wave optical cavity [3 [6] . In our 
setup, we can minimize parasitic diffusions processes [7] 
which may kick atoms out of the cavity and thus limit the 
lifetime T cav . The expansion of the atomic cloud, whose 
size quickly exceeds the diameter of the mirror, usually 
limits the number of bounces in the cavity [5J . As in [5J [B] , 
we circumvent this problem by using a "curved mirror" 
which refocuses periodically the condensate. This gives 
a very promising sensitivity for the proposed gravimeter, 
which increases as Tclv with the atom interrogation time 
T cav as in standard atom gravimeters £Q. 



I. DETERMINATION OF THE 
ACCELERATION OF GRAVITY: PRINCIPLE OF 
THE MEASUREMENT 

Following up our approach in [8], we present in this 
section a first simple description of the proposed matter- 



wave cavity and give a heuristic analysis of its perfor- 
mance as a gravimeter. 



A. Principle of the experiment 

The principle is to levitate a free falling atomic sample 
by providing a controllable acceleration mediated by 
a coherent atom-light interaction. Radiation pressure 
could provide levitation, but the resulting force is not 
precisely tunable if tied to incoherent spontaneous 
emission processes. A better choice to provide this 
acceleration is thus a series of vertical Raman pulses. 
Indeed these pulses impart coherently a very well defined 
momentum to a collection of atoms [9 . A sequence 
of Raman pulses of identical effective wave vector k 
interspaced with a duration T gives an acceleration to 
an atomic cloud which is monitored by the choice of T. 
Levitation occurs when the sequence of vertical Raman 
pulses compensates, on average, the action of gravity. 
This stabilization is obtained thanks to a fine-tuning of 
the period between two pulses: after a fixed time, one 
observes a resonance in the number of atoms kept in the 
cavity for the adequate period To. The atomic cloud is 
then well stabilized, and the average Raman acceleration 
equals the gravitational acceleration. Knowing the 
period To, one can infer the corresponding Raman 
acceleration and thus the gravity acceleration g. The 
ratio Hk/m can be simultaneously determined using the 
resonance condition of the Raman mirrors. 



B. Description of the cavity 

As displayed in Fig.l, the atomic sample, initially 
at rest in the lower state a, is dropped. After a free 
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fall of duration T/2, during which the sample acquires 
a momentum gT/2, we shine a first Raman tt pulse 
of 2 counterpropagating lasers with respective wave 



vectors kd c 



= ki and k up = k2 and respective 



frequencies oJdown = and ui up = u}%. This brings the 
atom from the internal state a to an internal state b 
with a momentum transfer 2hk = k% — k\. Then we 
shine a second Raman 7r pulse with kdown = k4 and 
k U p = k3, with respective frequencies LOdown = ^4 and 
uj up — u>3 (u>3 — u>4 ~ LU2 — u>±). This pulse brings the 
atomic internal state back to state a with an additional 
momentum transfer 2hk. If the two successive pulses 
are sufficiently close, this sequence acts as a single 
coherent "mirror pulse" which keeps the same internal 
state a and modifies the atomic momentum by Ahk. 
In particular, if the initial momentum is —2hk, the 
mirror simply inverts the velocity. This "mirror pulse" 
is velocity-sensitive: it reflects only the atoms whose 
vertical momentum belong to a tiny interval around 
a specific value pq. Thus, in order to bounce several 
times, the atoms should have the same momentum po 
immediately before each "mirror pulse" . This implies 
a resonance condition (JT|) on the period between two 
pulses. The adequate momentum p is set by the energy 
conservation during the pulse and fulfills the resonance 
condition rt2b. As depicted in Fig. [2j if the resonance 
condition dT| is satisfied, the sample will have a periodic 
trajectory in the momentum space. It is this periodicity 
of the atomic momentum for an adequate time-spacing 
T of two successive "mirror pulses" which yields the 
picture of a matter-wave cavity in momentum space. 
Fig. [3] represents an energy-momentum diagram of the 
atomic sample during a cavity cycle. 
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FIG. 2: Altitude and momentum of the atom sample as a 
function of time for the resonant period T — To and for dif- 
ferent initial velocities. The atom cloud always undergoes 
a periodic motion in momentum space even with a nonzero 
initial velocity. The figure illustrates three different cases of 
stable cavity in momentum space. 
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FIG. 1: Setup description. 



FIG. 3: Energy-momentum diagram of the condensate, zo 
represents the height from which atoms are dropped. 
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C. Resonance Conditions 

For the resonant period To, the momentum kick im- 
parted during each "mirror pulse" is equal to the mo- 
mentum acquired through the free fall between two such 
pulses: 



To 



Ahk 
mg 



(1) 



If the period T differs from To, the atomic cloud takes 
an acceleration a = Ahk/m (1/T — 1/T ). The average 
speed resulting from this acceleration drifts the momen- 
tum of the atoms from the optimum value po = — 2hku z 
satisfying the Bragg condition Q associated with elastic 
energy conservation. This drives the atoms progressively 
out of resonance and a part of the cloud will not be 
reflected by the "mirror pulses". When T differs from 
T), we thus observe a drift in position and momentum 
as well as a leakage of the condensate. It should be 
noted that if T = To, a non-zero initial velocity does 
not reduce fundamentally the number of bounces of an 
atomic sample (Fig. Indeed, only the periodicity of 
the momentum space trajectory matters, and provided 
that one adjusts the Raman detuning to account for 
the shift in the momentum po, the sample can still 
be reflected several times. The atomic sample merely 
drifts with a constant average vertical velocity, the only 
limitation being then the finite size of the experiment. 
Thanks to this flexibility in the initial velocity, the 
resonance observed is robust to an imperfect timing of 
the trap shutdown. 

Bragg resonance conditions express the energy conser- 
vation during the pulses: 



UJ 4 +UJ a + 



u) 1 + w b + Sac + 
(p + 4/ik) 2 



(p + 2hk) 2 
2mh 



2mh 



zmn 

uj 3 + uj b + 5 Ac + „_ t (2) 



2mh 

(p + 2hk) 2 



2mh 



p is the matter-wave average momentum immediately 
before the first Raman pulse, and Sac is the light shift. 
Thanks to this second set of conditions, which directly 
impacts the reflection coefficient, the "mirror pulses" 
act directly as the probe of the resonant time-spacing Tq 
expressed in ([TJ. When the atomic sample is dropped 
without initial speed, for a nearly resonant period 
T ~ T , the first pulse brings the sample at rest, so 
that both pulses play a symmetric role. If one does the 
replacements oj\ — + 073 and 0J2 — > W4, the mismatch in 
the two Bragg conditions is then equal in absolute value, 
yielding identical reflection coefficients for both pulses. 
One can then consider that the two Bragg conditions 
merge into a single one. We assume from now on that 
the atomic cloud has no initial velocity [i] , but the 



extension to the general case is straightforward. 
The two conditions (JTJ and ^ must be satisfied 
to ensure the resonance of the matter-wave cavity. 
Nonetheless, condition ([lj on the period is much more 
critical than the Bragg condition ([2]). Indeed, a slight 
shift in the period T from its optimum value To implies 
for the condensate an upward or downward accelera- 
tion. The increasing speed acquired by the atoms will 
generate, through Doppler shifting, a greater violation 
of the Bragg resonance condition and thus greater losses 
at each "mirror pulse" . Conversely, a mismatch in the 
detuning LU2 — u)\ with the adjustment T = To will only 
induce constant losses at each cycle. The observation of 
a resonance in the number of atoms, when one scans the 
period between two Raman pulses, is thus very sharp 
and robust to an imperfect adjustment of the Raman 
detuning. Consistently, we choose to determine the 
acceleration of gravity g through condition 0. Fig. [4] 
sketches the number of atoms in the cavity as a function 
of the period T after different numbers of cycles with a 
detuning matching perfectly condition ([2|. We observe 
that the resonance in T becomes sharper as the number 
of cycles increases. 
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FIG. 4: Number of atoms in the cavity after 1,10,50 cycles as 
a function of the ratio T/T (To ~ 3.8 ms). We took f2 = 
2-7T x 5 x 10 3 Hz. The number of atoms initially present in the 
cavity was fixed to be N = 10 6 . 



D. Expected Sensitivity 

Let us derive the resonance figure associated with the 
period T. We characterize the mismatch in the Bragg 
condition by an "off-Braggness" parameter y(p) function 



[i] Bose Einstein condensates can be brought at rest very accurately 



and are thus well suited for our system. 
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of the average momentum of the atoms [TO] : 



y(p) 



2n 



[hk 2 /2m+ 2p • k/m — (wi — W2 — u ba — Sac)] 

where f^o the effective Rabi frequency of the Raman pulse 
and p the momentum immediately before the first pulse.. 
The Raman detuning is adjusted to be resonant if p = 
PqU z = —2hku Zl so we adjust the detuning to match 
y{-2hku z ) = 0: 

= (p z + 2hk)k = (p z + mgT /2)k 
2mf?o 2mf2o 

In the remainder of this section, we focus on the vertical 
component of momentum which we note p to alleviate 
the notations. Because of the mismatch in the Bragg 
condition, only a fraction p(p) of the atomic cloud will 
then be transferred during the first Raman pulse: 



Pip) = 



sin 



2 / 7T 



Wl+yip) 2 ) 



1 + y{pY 



(5) 



For the second Raman pulse, the mismatch is equal in 
value and opposite in sign, so the same fraction of atoms 
undergoes the second transition. The reflection coeffi- 
cient of the "mirror pulse" is then simply the product of 
these values: 



R( P ) 



sin 



(l + y(p) 2 ) 2 



(6) 



The part of the cloud which is not reflected will simply 
go on a free fall and have an off-Braggness parameter of 
y{p — mgT) for the next "mirror pulse" . We will adopt 
experimental parameters such that yimgTo) 3> 1, so that 
non-reflected atoms arc insensitive to subsequent Raman 
transitions and can be considered as expelled from the 
cavity. 

The average momentum acquired by the atomic cloud 
results from a competition between the gravitational ac- 
celeration and the kicks of the "mirror pulses" . Given a 
period T for the sequence, the average vertical momen- 
tum of the cloud immediately before the n-th "mirror 
pulse" is simply: 



Pn 



-mgT j2 + (n - 1) x mgiT - T) 



(7) 



The remaining fraction of the cloud after n cycles is thus: 

R{T) = Rip x )...R(p n ) (8) 

We expand this expression in Appendix [X] for nearly res- 
onant pulses. We have assumed that a relative variation 
e of the condensate population can be tracked experimen- 
tally. This computation shows that the error in the de- 
termination of the gravity acceleration can be less than: 




JW-log(l-e)' 



3/2 



(9) 



where v r = hk/m is a recoil velocity which can be mea- 
sured independently or, as stated before, directly from 
resonance condition ^ [TT]. This velocity has been de- 
termined with an accuracy as good as a few 1CP 9 for 
Cs [T2J [T3] and Rb atoms [J3] , ultimately limits the per- 
formance of our gravimeter. Fig. [5] displays the relative 
error in the determination of the acceleration of gravity 
obtained from a numerical simulation with the reflection 
coefficient (pj. It shows a very good agreement with the 
analytic expression (|9| after about 10 cycles. 
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FIG. 5: Error on Ag/g as a function of the number of cycles, 
with Raman pulses of Rabi pulsation fio = 2-7T x 5 x W 3 Hz 
and a detection threshold of e = 10 -3 . The lower curve is the 
actual sensitivity based on a simulation with a condensate of 
infinitely narrow momentum distribution. The upper curve is 
the analytic formula jjjj. 

Formula |9| thus yields a sensitivity which scales as 
n 3/2 _ iT cav /T ) 3 / 2 , where T cav is the interrogation 
time of the atoms and To the duration of a cycle at 
resonance. We thus obtain the expected improvement of 
the sensitivity with the atom interrogation time. This 
is normal since the accuracy of the measurement also 
increases with the selectivity of the Raman pulses. The 
mirror pulses should thus be as long as possible to act 
as efficient atom velocity probes. Ideally, each pulse 
should last half of the optimal period To. The atom 
sample would then fall in a continuous light field. This 
setup ressembles that of Clade et al. [3], except that 
here the atomic sample is interacting with a travelling 
wave, and that in addition to condition (JT|) we have a 
double Raman condition d2j. 

In the previous description, one could in principle 
maintain the sample for an arbitrary long time inside 
the cavity provided that the two resonance conditions 
are fulfilled, yielding an extremely accurate measure- 
ment of the gravitational acceleration. In fact, even at 
resonance, systematic losses occur that limit the number 
of condensate cycles. Practically, one can hope to keep a 
significant number of atoms in the cloud up to a certain 
number of n cycles, which reflects the sample lifetime 
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T 



nTo at resonance. These losses at resonance have 



several origins. 

First, they can result from imperfect recoil transfers 
due to residual fluctuations in the Raman lasers inten- 
sities. Raman impulsions can be very effective though, 
since these transfers have been realized with an efficiency 
as high as 99.5%. 

Second, parasitic diffusions may eject atoms from the 
cavity. These processes, such as absorption followed by 
spontaneous emission, can be made arbitrarily small by 
using far-detuned Raman pulses. This limitation is thus 
essentially of technical nature. 

Third, the cloud expansion can drive the atomic sam- 
ple out of the Raman lasers. Indeed, if the cloud is not 
refocused, its transverse size exceeds quickly the diameter 
of the Raman lasers. The atoms would then be limited to 
a few cycles in the cavity. It is thus essential to involve a 
mechanism that stabilizes transvcrsally the atomic cloud. 
We expose in Section |III| two ways to focus the atomic 
sample. 



II. ABCD ANALYSIS OF THE MATTER- WAVE 
CAVITY 

In the following, we will calculate explicitly the evolu- 
tion of an atom sample in our gravitational cavity using 
the ABCD matrix formalism [T5J [Tj)J [T7] . In this section, 
we assume that the atom density after the initial free 
fall is sufficiently low to make the effect of interactions 
negligible during the subsequent bounces. 



A. Description of the atomic sample 

We will restrict ourselves to a Bose-Einstein conden- 
sate which, in the Hartree-Fock approximation, can be 
described by a macroscopic wavefunction ^(r,t). After 
the initial free-fall, the evolution follows the linear partial 
differential equation: 



mMr,t)) = {p 2 /2m + H G }\*(r,t)) 



(10) 



where Hq is the Hamiltonian associated with gravito- 
inertial effects. From the linearity of this equation, we 
can deduce the time evolution of any arbitrary wave func- 
tion from the propagation of a complete set of functions. 
As explained in Appendix [Cj the propagation of such 
a basis, the Hermite-Gauss modes Hi mn (r), can be ex- 
tracted from the propagation of a generating Gaussian 
wave function |17j : 



*7?(r) 



exp[^(r-r )Y Xo\r-r ) 



+ i(r - r„) • (Po - 2hxl 1 lZ) + \~^X^X^] (11) 

The matrices X ,Y and the vectors r ,Po correspond 
respectively to the position and momentum widths, to 



the average position and to the average momentum of 
the wave function. We can then restrict our deriva- 
tion without any loss of generality to the propagation 
of (111. Since the Hamiltonian Hq can be considered 
with a very good approximation to be quadratic in po- 
sition and momentum, the wave-packet (111 follows the 



ABCD law for atom optics [IB]. In order to alleviate the 
notations, we shall omit to mention the index ~a in the 
subsequent computations and denote the corresponding 

State |*ro,po,Xo,Yo)- 



B. Initial Expansion 



Before the trap shutdown, the condensate evolves un- 



We 



dcr the Hamiltonian H = p /2m + Hq + -ffxrap- 
remove the gravitational term Hq thanks to a unitary 
transform UQ(t,to): 



m)) = u G (t,t )\m) 



(12) 



UQ{t, to) represents the evolution of a quantum state un- 
der a gravitational field. Performing this unitary trans- 
form is equivalent to study the condensate in the non- 
inertial free falling frame. The state |3>(i)) then evolves 
under the Hamiltonian Hp ree p a u = p 2 /2rn+_ffrp ra p. The 
condensate is taken to be initially in the strong cou- 
pling regime, so that the corresponding wave function 
<&(r, t) follows the scaling laws established by Castin and 
Dum 18J for the Thomas-Fermi expansion. The initial 
wave function <&(r,i ) corresponds to the Thomas- Fermi 
profile. We represent concisely its evolution by the uni- 
tary transform UTFE(t,to): 



Mt)) = U T FE(t,t Q )\$(t )) 



(13) 



It is indeed not useful at this point to explicit this trans- 
form, whose expression is given in Appendix[B] We trans- 
form back to the laboratory frame at the time t\ when 
we start to shine the first "mirror pulse": 

|tf(tl)) = M*l,*0)tfTF B (tl,t0)|$(*0)) (14) 

The resulting quantum state ^(fi)) will be taken as a 
starting point for the subsequent oscillation of the con- 
densate in the cavity. Its evolution is conveniently ob- 
tained by decomposition on a suitable basis of Hermite- 
Gauss modes Hi mn (r), as explained in the preceding 
paragraph. The initial free fall simply determines the 
initial coefficients of the projection: 



*(r,ti) 



l,m,n 



(tx)Hi mn {r) 



(15) 



C. Propagation of a Gaussian wave-packet in the 
diluted regime 

The evolution of free falling atoms in a Raman pulse 
is non-trivial since the gravitational acceleration makes 
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the detuning (|3| time-dependent: 



A(f) = CJl — UJ2 — LUba 



2(p — mg t u z ) • k/m 
-2/ik 2 /m - Sac 



(16) 



The behavior of a two-level atom falling into a laser wave 
has been solved exactly [TO]. Gravitation alters signifi- 
cantly the two-level atom state trajectory on the Bloch 
sphere when the pulse duration t becomes of the order 
of: 



1 



v4W 



10" 4 s 



(17) 



Indeed, for this duration the off-Braggness parameter ^ 
changes significantly during the pulse. As seen in the pre- 
vious section, in order to probe effectively the resonance 
condition ([TJ, the Raman 7r-pulses need to be velocity- 
selective. This leads us to consider pulse durations on 
the order of the millisecond, typically longer than r g . 
It is then necessary to compensate the time-dependent 
term induced by the acceleration of gravity in the de- 
tuning by an opposite frequency ramp chirping the pulse. 

The simultaneous effects of gravito-inertial and elec- 
tromagnetic fields can be decoupled thanks to an ef- 
fective propagation scheme developped by Antoine and 
Borde [16j EZ] ■ It accounts for the electromagnetic inter- 
action through an instantaneous diffusion matrix S and 
for gravito-inertial effects through a unitary transform 
C/i(T,0): 



|tf(T)} = tfi(T,0)S|tf(0)) 



(18) 



Following this propagation method, it is sufficient to ap- 
ply an effective instantaneous diffusion matrix for each 
"mirror pulse" and evolve the state between the pulse 
centers as if there was no electromagnetic field. 



D. Action of the effective instantaneous interaction 
matrix S 

We study in this paragraph the interaction of the con- 
densate with a quasi-plane electromagnetic wave, for 
which the instantaneous diffusion matrix S is known. 
This matrix is operator- valued, but momentum opera- 
tors can be taken as complex-numbers since the consid- 
ered wave function is a narrow momentum wave-packet 
centered around a nearly resonant momentum po (i.e. 
such that y(po) "C 1). In other circumstances, this inter- 
action can give rise to fine structuring effects such as the 
splitting of the initial wave into several packets following 
different trajectories (Borrmann effect) [JO]. Following 
the approach of the paragraph |II A| we consider a Gaus- 
sian matter- wave: 



l*o) 



la,* 



We study the interaction of this atomic wave with a "mir- 
ror pulse" involving two linearly polarized running laser 
waves: 



E = E(x, y) cos (kz - uit + 
+E(x, y) cos (kz + uj 2 t + $ 2 ) 



(20) 



With respect to the population transfer, electromagnetic 
fields may be treated as plane waves in the vicinity of 
the beam waist. The effective diffusion matrix 5(k,p ) 
associated to a Raman pulse effective wave vector k and 
applied to a wave-packet of central momentum po then 
yields [10] : 



5(k, Po ) 



Sbb 
Sab 



Sba 
Saa 



(21) 



Saa — Sfrb — exp 



-t(n AC (v)+n AC (v))- exp(i<5 12 r) 



cos 



(fl T y/1 + y- 



+ 1- 



S a b = S ba = i exp 



— ? - in (^otV 1 + y 2 ) 

i(Q e AC (r) + Q b ac (t)) J] exp (-i6 12 r) 
x sin (q t y/1 + y 2 ) /y/l + y 2 



2/(Po,k) 



-(5i 2 (p ,k) 
2O 



1 



[lui - uj 2 - u a b - k • Po/m 



2fl 1 

-fik 2 /2m-^ c (r)-0» c (r))] 

(22) 



u>i , Ll>2 are respectively the pulsations of the lasers prop- 
agating upward and downward, r the duration of the 
pulse, £l AC (r) and fl AC (r) are the AC Stark shifts of 
the associated levels. It is worth commenting the posi- 
tion dependence of those terms, which intervene in two 
different places in the S matrix. In the off-Braggness 
parameter y, the term SAc( r ) induces an intensity mod- 
ulation, while in the complex exponential, the term 
f2o(r) = Sl AC (r)+Sl AC (r) changes the atomic wave-front. 

After each "mirror pulse" , the part of the condensate 
which does not receive the double momentum transfer 
will fall out of the trap if y(po + mgT) 3> 1. We thus 
project out those states and focus on the non diagonal 
terms of the diffusion matrix: 

|*(2r)) = <a|5(-k,p o +?ik)|6)(6|5(k,po)|*(0)) (23) 

The state after the mirror pulse is thus: 

|*(2t)) = /o(r)|a,* ro , Po+2ftk ,x ,y ) 



(fWi + y(po) 2 ) 2 



(24) 



The amplitude factor p(r) reflects both the loss of non- 
reflected atoms and the change in the atomic beam wave- 
front. Expanding the generating wave function (11) into 



powers of a, one shows that the effect of the interaction 



(19) matrix S is the same on each mode of the expansion ( 15 1 
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E. Gravito-inertial effects 



A. Matter-wave focusing with phase imprinting 



The unitary transform U±(T, 0) represents the gravito- 
inertial effects. We refer the interested reader to [17] 
for a thorough derivation of this operator. We remind 
here the main result necessary for our computation. This 
operator maps a state defined by a Gaussian of parame- 
ters Xi ,Yi,ri,pi in position representation onto an other 
Gaussian state in position representation whose parame- 
ters X2, Y2, T2,P2 depend linearly on the former according 
to: 



cosh^T] 7-V2 sinh^T] \ / Xi 
7 1 / 2 sinh^T] cosh^T] J 



X 2 

Y 2 J V 7 1/2 sinhk/7Tl coshf^Tl J \ Y x 

(25) 

The coefficient 7 reflects the interaction effects through 
an effective potential quadratic in position and assumed 
to be constant in time. There is the same matrix rela- 
tion between the initial and final position and momentum 
centers ri,pi,r2,P2 of the wave-packets, with an addi- 
tional function £ which reflects the constant part of the 
gravity field. The transform £/i(T, 0) also introduces an 
additional phase factor given by the classical action: 



J7i(T,0)|a,* 



ri,pi,Xi,Yi, 



= e iS C i(T,0) 



a,* 



r2.P2,^2,V 2 / 



Since this phase factor does not play any role in the fol- 
lowing computations, we do not give its expression here, 
but it can be found in reference |16j . Expanding a generic 



Gaussian such as ( 1 1 1 shows that the propagation of any 



Hcrmite mode of the expansion ( 15 1 is identical in the 
gravity field: the Gaussian parameters X, Y involved in 
each mode are transformed identically. 



F. Conclusion: cycle evolution of the matter wave 

In our approach, the effect of the interactions has been 
neglected after the first bounce and the propagation of 
the diluted atomic sample in the cavity is essentially 
mode-independent. Nonetheless, in experiments where 
atomic samples of higher density are bouncing on electro- 
magnetic mirrors jSJ HHJ |5TJ , interactions do change 
the shape of the cloud during the propagation. As we 
shall see in the next section, interactions impact the 
transverse velocity distribution in a way that can lead 
to a reduced stability of the cavity. In the following, 
we will review possible focusing techniques to solve this 
problem. 



III. MATTER- WAVE FOCUSING 

We investigate in this section two possible curved mir- 
rors. We first review a focusing technique based on 
the phase imprinting through a position-dependent Stark 
shift [23] - Afterwards, we introduce an original focusing 
mechanism based on a laser wave-front curvature trans- 
fer. 



This method has the advantage of leading to tractable 
equations. It relics on a position dependent Stark shift 
provided by quasi-plane waves with a smooth intensity 
profile: 



E(x,y) = E Q [l- 



y 



(27) 



Unfortunately, this Stark shift implies a position- 
dependence due in the population transfer. This results 
in a loss of atoms which makes this focusing method 
hardly compatible with the extreme cavity stability re- 
quired by this experiment. Nonetheless, it is interesting 
to demonstrate the effect on the wave curvature induced 
by this position dependent light shift. In this perspec- 
tive, we neglect the position dependence in the popula- 
tion transfer but not in the phase of the diffracted matter 

wave. Indeed this shift Q A c(r) = V-° AC (l - 2^j^) im- 
prints a quadratic phase to the matter wave: 



|*(2r)) =piexp[- 



(x 2 +y 2 )}\a,^ To 

,Po+2fik,X ,io/ 

(28) 



(26) 

The position-independent phase shift is hidden in the co- 



efficient pi. Using expression (111 for the wave function 



v ^ro,po+2fik, J f ,y 0J one can recast the last equation into: 



|*(2t)) =pi|a,* ro , Po+2 nk,x 1 ,r 1 ) 



with: 



SMV\ 



X 
Y 



(29) 



(30) 



l AC T p, 

....2 r± 

with Pj_ projection matrix on the transverse directions 

/l 00 
P x = 1 Q 
\0 

The AC shift factor thus changes the Gaussian parame- 
ters of the matter wave just like a thin lens of focal / in 
classical optics, where the transform law yields: 



Xi 

Yj 



1 



1 



X 
Y 



(31) 



One could thus define the focal length of an atom optic 
device as the / parameter entering the ABCD transform 



(30 1. Precisely, a phase factor exp [— ia(x 2 



y 2 + z 2 )] 

on a Gaussian atomic wave changes the Gaussian pa- 
rameters of the matter wave according to the ABCD law 
of a thin lens of focal: 



/ = 



2ha 



(32) 



The AC dependent Stark shift thus plays for the atomic 



2 /8hn° A 



AC' 



beam the role of a thin lens of focal / = 
Let us point our that this focusing occurs in the time 
domain so that the "focal length" is indeed a duration. 
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B. Matter-wave focusing with spherical light waves 

As mentioned in the last paragraph, the impossibility 
to maintain a perfect population transfer on the whole 
wave-front while focusing with a light shift effect makes 
this technique inadequate for the proposed gravimeter. 
We investigate here a different method which does not 
have this major drawback. Instead of shaping the atomic 
wave-front thanks to an indirect light-shift effect, a better 
way to proceed is indeed to have the matter wave interact 
with a light wave of suitable wave- front, like on Figure [6] 
Following up this intuitive picture, we propose an alter- 



Matter wave phase front Raman wave phase front 




Refocalized matter wave 
phase front 




FIG. 6: Interaction between a spherical matter wave and a 
laser. The laser spherical wave-front refocuses the matter 
wave. 

native matter- wave focusing scheme, fully original to our 
knowledge, based on the matter-wave interaction with 
electromagnetic fields of spherical wave-front. To show 
that the focusing is effective, we compute the transition 
amplitude of a matter wave interacting with Gaussian 
Raman waves to first order in the electromagnetic field. 
A similar computation has been previously performed by 
Borde in the context of atomic beamsplitters (24j . 

Let us consider for the "mirror pulse" two counterprop- 
agating matched Gaussian beams: 



photon into mode Ei . We have used the confocal param- 
eter of the light beam b = k^w 2 as well as the complex 
Lorcntzian function L + 25 1: 



Ei(x,y,z,t) = |E#(r-i 
E 2 (x,y,z,t) = iJ7 "(r-r 



\ p ik a (z-z w 



E(t)e 
)E{t)e 



L + (z) 



c.c 



1 



1 — 2iz/(fc Wo) 



^=^= exp[i arctan f— 1] 

i I z 2 ^ ZO/ 

1 4b 2 



The term arctan (^) is known as the Gouy phase. The 
matching of the two laser beams is reflected in the 
relation between their transverse structures U~(r) — 
U + *(r): at each point their curvature is identical. The 
Raman diffusion associated with these fields yields an 
effective interaction Hamiltonian whose matrix elements 



Vba(r,t) = -fifi U+(r - r w )f/ -*(r - r w )e 



i2ko(z— z w ) 



xe 



-i(ui 2 i+r(t-t r ))t+iip 



C.C. 



F(t) 
(35) 



The term r(t — t r ) accounts for the frequency ramp start- 
ing at time t ri UJ21 — — ^i is the Raman detuning 
and F(t) = \E(t)\ 2 the time envelope of the pulse. The 
computation of the transition amplitude is somewhat in- 
volved and deferred to Appendix |D| We obtain: 

xL+(z - z C0 (t))U+ 2 (r - rco(t) - — (t - «o) 
\ m 



x(b,v\U (t,t ) 



d 3 p 

(2ttH) 3 / 2 



F(w B (p„ko))<o,p|*(to)) \b,P 



(36) 



Uo(t,to) is the evolution operator in the gravitational 
field, r co (t) = r w +^(t-t ) + ig(t-to) 2 and iv B ( Pz , k ) 
a frequency which reflects the Bragg resonance condition: 



Ws(p z ,k ) = -[Wba 



2m 



W21 



:(*r 



-to)) 
(37) 

^^"^(r) corresponds to a Gaussian mode of confocal 
parameter b — fcou>Q and waist wq/\/2. The first-order 



term ( D32 1 is the leading contribution to the outgoing 
excited matter wave. The filtering of the pulse acts as 
expected through the Fourier envelope F (ujb(p z , k )), 
significant only for a small velocity class which can 
be tuned by the starting time t r of the velocity ramp. 
The operator Uo(t,to) reflects the propagation in the 
gravitational field. The factor L + (z — Zc(t)) barely 
affects the longitudinal shape of the atomic wave, 
without contributing to the average vertical momentum. 
As discussed in Section [TJ the cavity lifetime of the 
atomic sample does not depend on its longitudinal pro- 
file. Therefore we do not need to worry about this factor. 



(33) 



What matters is the transverse structure of this outgo- 
where fco = (k2 — ki)/2. The detuning 102 — ^1 is adjusted ing wave, which corresponds to a focusing matter wave, 
so that the relevant Raman process be the absorption of As suggested in Figure [6] the curvature of the Gaussian 
a photon from mode E2 followed by the emission of a Raman wave has been transmitted from the laser wave 
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to the atomic wave through the term Uq [r — rco(^) 



2Mc 



^(t — t )]. This term induces a quadratic dependence 
of the phase on spatial coordinates: 



U+'(r) = exp[- 



i4k a z 



k$w$ + 4z 2 



(x 2 + y 2 )} 



(38) 



Following the approach of the precedent paragraph, and 
the relation ( 32 ) , this can be interpreted as a thin lens 
effect. To express the corresponding focal, we introduce 
the vector f (t) = r — r C o{t) — — *o)- For an atomic 

wave centered around z = z a at the time i, the interaction 
with the light field plays the role of a thin lens of focal 
f(z ai z w ,t) to first order in the electromagnetic field: 



m kZw. 



8k z(t) 



(39) 



where the parameter z(t) is: 



POz 



2hk 



(t-to)- 7;9{t-hf 



(40) 

This focal can thus be controlled by the relative position 
of the laser waist and atomic wave-packet center. This 
first-order computation shows how the laser beam cur- 
vature is transferred to the matter beam wave-front and 
suggests that it is possible to focus matter waves through 
Raman pulses with a spherical wave-front in a control- 
lable manner. 



C. Transverse stability of the cavity 

The main threat to the cavity stability is indeed the 
interatomic interactions which will push away the atoms 
of the sample from the central axis of the cavity. It is 
indeed possible to approximate these interactions with 
an effective lens. A detailed ABCD matrix analysis of 
interaction effects is given in |26j . Here we just consider 
that interactions induce an effective quadratic potential 
represented by the diagonal matrix 7^. 

In the precedent paragraphs, the focusing obtained is 
only effective for the transverse directions. This is suffi- 
cient, since the longitudinal spread of the atomic sample 
does not drive it out of the laser beams. The transverse 
stability of the cavity is entirely reflected in the temporal 
evolution of the Gaussian parameters X(t),Y(t). The fi- 
nal parameters are related to the initial parameters by 
the ABCD matrix: 



cosher] 



i 



: sinh[^77T] 



yrTsinh^ 
cosh[, 



%T] 7 .r 1/2 sinner] \ 

/7IT] cosher] ) 

(41) 

As in |26) , one can use this input-output relation to model 
the stability of the matter-wave cavity. For the diluted 
matter waves involved in this system, a slight curvature 
in the mirror is sufficient to stabilize transversally the 
atomic beam. 



IV. CONCLUSION 

We would like to point out the physical insight 
provided by the picture of a cavity in momentum space. 
Such a cavity is only possible in atom optics since 
photons, whose velocity is fixed to c, cannot be accel- 
erated. In our system, the corresponding momentum 
wave-packet oscillates between two well-defined values 
(Fig. |2j, with a resonance observed for the adequate 
time-spacing of the mirrors. We can push further the 
analogy with an optical cavity. The force mg is the 
speed of the field in the momentum space. 4hk is the 
momentum analog for the cavity length. The cycle 
period Tq of the wave propagating in momentum space 
is orders of magnitude longer than the usual cycle time 
of a pulse in an optical cavity. 

Let us look again at the resonance condition with 
this picture in mind. This relation corresponds in mo- 
mentum space to: 



T = L/c 



(42) 



with the replacements mg — > c (propagation velocity in 
momentum space) and Ahk — > L (distance in momentum 
space) . The usual resonance condition for an optical cav- 
ity yields: 



T = nL/2c neN 



(43) 



The difference can be explained by the fact that, unlike 
in an optical cavity where the light goes back and forth 
between the mirrors, the "way back" in momentum 
space from — 2hk to 2Kk has to be provided by the 
light pulse. This explains why the factor 2 is absent 
in the denominator of ( 42 ) . The integer n is absent 



in the resonance condition ( 42 1 because we considered 



only two-photon Raman pulses for the optical mirrors. 
Indeed, for each period T n — hTq, the cavity becomes 
resonant with mirror pulses based on n-photons pro- 
cesses. In order to levitate, the atomic sample needs 
to receive from the light pulse an adequate momentum 
transfer of Anfik.. This momentum fixes the number of 
photons exchanged for each possible resonant period T n . 

In the proposed gravimeter, the momentum cavity 
is loaded with short single atomic pulses well-localized 
in momentum space, since the instantaneous velocity 
distribution is sharply peaked at any time. This is the 
analog of a femtosecond pulse propagating in an optical 
cavity. It would be interesting, however, to load the 
cavity with a continuous flow of free falling atoms coming 
from a continuous atom laser. At a fixed momentum, 
the contributions from different times would sum-up and 
interfere, exactly like in a Perot-Fabry interferometer. 
This system would then constitute to our knowledge the 
first example of a momentum space cavity continuously 
loaded with a matter-wave beam. 
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We have studied the levitation of an atomic sample 
by periodic double Raman pulses. In our system, the 
matter wave is trapped in an immaterial cavity of peri- 
odic optical mirrors. For the adequate time interspace 
between two pulses, the atomic sample is stabilized and 
levitate for a long time. Thanks to the sensitivity of 
the stabilization to this period, one obtains an accurate 
determination of the gravitational acceleration. 



has y(pi), y(p n ) 1. The expression for the reflection 
coefficients simplifies to: 



R(Pi) = 



1 



(i + y 2 fe)) 2 



+0(y\ Pi )) = l-2 2/ 2 (p l )+0(y 4 ( Pl )) 

(A3) 

The fraction of atoms kept in the cloud can then be ex- 
pressed as: 



In our approach, the system could be loaded with 
any atomic sample describable by a macroscopic wave 
function. It is indeed not necessary to impose an initial 
small velocity dispersion, since the first mirror pulse 
will serve as a filter for a narrow velocity-class, while 
the next pulses will serve as a probe. Many aspects 
developed in this paper are still valid for a thermal 
cloud. Nonetheless, Bose Einstein condensates are 
ideally suited for this trap since matter-wave focusing 
is more efficient with a single mode coherent source. 
In this paper we have considered only 7r-pulses for the 
atom- light interactions. In fact, one could consider 
other schemes, for example one could split each 7r-pulse 
in two copropagating 7r/2-pulses separated by a dark 
space resulting in a sequence of Ramsey-Borde inter- 
ferometers [27] . Since the sensitivity to gravitation is 
proportional to the area covered in space-time by the 
interferometer, the optimal situation is obtained when 
the copropagating 7r/2-pulses are separated by To/2. An 
experimental realization of this proposal is planned with 
the support of the Institut Francilien de Recherche en 
Atomes Froids(IFRAF). 
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APPENDIX A: COMPUTATION OF THE 
GRAVIMETER SENSITIVITY 

After n nearly- resonant cycles, the fraction of the cloud 
preserved becomes: 



R(T) = R( Pl )...R(p n ) with R(p) 



(i + y(p) 2 ) 2 
(Al) 



with expression ^ for the momenta: 

Pn = - mg T/2 + (n - 1) x mg(T -T) (A2) 

At resonance T — To, one would have pi = ■■ = p n = 
—2hk and y(p\) = .. = y{p n ) = 0. For a big number of 
cycles n, in the vicinity of the resonance T ~ Tq one still 



1 It n 

log ^ = - E lo s ( fl (p*)) - E lo s C 1 + 2 ^ 2 ^)) 

^ ' i=l i=l 

(A4) 

We insert expression ([7]) for the momentum in Q to 
derive an expression for the off-Braggness parameter: 



, \ (j(T - T )k . 

y(Pi) ^ — o x 1 

"0 



(A5) 



The reflection coefficient becomes: 



Ml 



This sum may be approximated by an integral because 
i > 1: 



1 



log 7 ^ v ~ / dx log ( 1 + 2x 2 b ^ ^; u/ - ) (A7) 



!ff 2 (T-T ) 2 fc 2 



R(T) ' ji 

This integral can be performed analytically: 

dx\og(l + ax 2 ) = ~2x+2 ^ ^ + xlog(l + aa; 2 

V a 

(A8) 

We set a = 2 5 2 (T - T ) 2 fc 2 /fig, which verifies y/an ~ 
y(p n ) <C 1. We can then use expression (A8) in (AT) and 
Taylor expand the right hand side: 

log -J— ~ -2n + 2[n - \arV i ] +nan 2 = \an :i (A9) 
Kyi ) o o 

We have omitted the small term coming from the lower 
bound of the integral. We finally obtain: 



log 



1 2 5 2 fc 2 (T-T ) 2 n 3 



R(T) 



3ft 2 



(A10) 



Let the condensate perform n bounces for a range of 
values of T close to the expected value To, and de- 
tect the number of atoms in the cloud afterwards. If 
a relative variation e can be tracked experimentally, we 
can bound the period To between Ti and Ti such that 
R(Ti) = R(T2) = 1 — e. According to our previous com- 
putation: 



|3 _ r>|i|f| ... J^f^l^ (A11 > 
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We infer the gravitational acceleration from the period 
To thanks to relation ([I]) , so that their relative errors are 
related by: 



l^l < if l 

9 T 



AVr 



(A12) 



with the recoil velocity v r = hk/m. This gives the follow- 
ing upper bound for the relative error on the gravitational 
acceleration g: 



q " V 8hk 2 



'fW-log(l-e)' 



,3/2 



(A13) 



APPENDIX B: THOMAS-FERMI EXPANSION 

The evolution of a condensate initially in the strong 
coupling regime yields [T8] : 



U T FE(t,t )Mto)) = 



(Bl) 



with: 

and$(r,i )^ ( ] 4) 1/2 (l-ci 



VAi(t)A 2 (t)A 3 (t) W JW ' 



(^ 2 +y 2 ) ; .2 Z ^ 1/2 

_R2 Z 2 



For a cigar-shaped condensate, the frequency ratio e = 
^p- is small and we may keep track of the radial expansion 
only: 



A,(t) = l X±(t) = y/l + u** 
/?(<) = arctan (l + w 2 * 2 ) 



(B2) 
(B3) 



APPENDIX C: PROPAGATION OF A WAVE 
FUNCTION: THE METHOD OF THE 
GENERATING FUNCTION 

Let us assume that we know the solution of a linear 
PDE for a family of initial conditions indexed by ~a: 



d t f(r,t,~a) = L(/)(r,i,"a) 
: \73efe) exp[ ^' 



/(r, 0, 7?) = - 7T h^ exp[ff (r - r ) - r ) 



From the linearity of L, the propagation of an Hermite- 
Gauss mode Hi mn , i.e. the solution at future times of 
the partial differential equation with initial condition: 



d t gimn{*,t) = L(gi mn )(r,t) 
gi mn (r,0) = Hi mrl (r) 



(C3) 



can be inferred from the coefficient of a\u™al{ in the a 
expansion of /(r, t, ~a), after a change of variables in the 
argument of the Hermite polynomial. The propagation 
of any arbitrary wave function ^>(r,t) then follows by 
linearity from the computation of the initial projection 
on the Hermite-Gauss basis: 

*(r,t) = ^ c lmn9lmn(*,t) 
with c lmn = J d 3 rH; mn (r)y(r,0) (C4) 



APPENDIX D: COMPUTATION OF THE 
FIRST-ORDER TRANSITION AMPLITUDE 
WITH SPHERICAL WAVES 

The state vector evolves under the Hamiltonian H = 
Hq + He + V, where Hq accounts for the internal atomic 
degrees of freedom, He = p 2 /2m + mgz for the external 
particle motion and V for the light-field. In this appendix 
we compute the transition amplitude to first order in 
V. To perform this computation, we consider the state 
vector |^(t)) in the interaction picture: 



(Dl) 



Uo(t,ti) is the free evolution operator in the absence of 
light field between times t\ and t: 

U (t, ti) = exp (-iH E (t ~ h)/h) exp(-vff (t - h)/h) 

(D2) 

The light field is turned on at time to- The first-order 
term of the Dyson series associated with the potential V 
is: 

\* {1) (t)) = \ f dt'V(t')\¥°\t )) (D3) 
where V(t') is the potential in the interaction picture: 



-f(r-r o )-(po-2^~- 1 ^) + i^X - 1 X *^] = * 1? (r) V(t) = U^(t, i x ) (Ha(r op , t) ® \b) (a\) U (t, t x ) + h.c. 



(CI) 



where L is a linear differential operator in the first two 
variables of the function /. Hermite modes can be defined 
through an analytic expansion of the exponential ^!-g(r): 



exp[rMir + vM2~3 + a M3 a] 
= erMir Ei, m ,n i l+m+n a[a^a^H lmn (M 2 r,-^M 3 ) 



(C2) 



^V ba (R op (t,t 1 ),t)<»\b)(a\e^ t - t ^ + h.c. 

(D4) 

R-op(^^i) is the position operator in the interaction pic- 
ture, given by integration of the Heisenberg equation of 
motion |17j : 

Rop(Mi) = ^(MiKpt/otMi) 
= A(t, h)r op + B{t, h)p ov + ^(i, h) (D5) 
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The parameter t\, associated with a choice of repre- To compute the Fourier components of V, we first use 
sentation for the interaction picture, can be chosen as the BCH relation: 
t\ = t . We need only consider the term V ba of the inter- 
action potential, for which we adopt the usual rotating 



wave approximation. The first-order transition ampli- 
tude &W(r, t) is given by the relation: 



x e 2 



ik-(^(t',to)r op +B(t',to)Pop+£(*',to)) _ e ik-A(t',t )r O p e ik-B(t',t )po 
|[A(t',to)k-r p,B(t',to)k-p op ]+iC(t',to)) 



6«(r,t) 



i (b,r\U {t,t ) f* o dt'V ba (R OI> (t',to),t') <g> |6)(a| 



xe— ( t - t »)[/ - 1 (< ,<o)l*(io)) (D^ 



In order to understand how the light wave structures the 
atomic wave-packet, we introduce the matrix elements of 
V between plane atomic waves: 



b(r,t) = -(b,r\U (t,t ) 
in 



df 



, dpdp' 



|b )P ') e ^(*'-*o) 



to ( 2 ^) 3 

x (p'|Ha (Ro P (t', to), |p)e^(*'- t0 ) (a, p|*(t„)) 



(D7) 

We introduce the Fourier transform of the interaction 
potential: 

V ba (r,t) = -nnF{t)e^ {t ' u)2 ^ iuJ2i{t - to)+lipo 

d 3 k 



(2tt)3/: 



: W(k)e 



ik-(r— r w ) 



(D8) 



he last commutator is responsible for the recoil term, 
and can be written: 



^[A{t',t )k-r op ,B(t',t )k-p op ] = l —kA(t',t )B(t',t )k 

( D1 °) 

The matrix element of the interaction potential contains 
the following term: 



^p'|^/(J t ) e ik-(A(t',to)r op +B(t',t )p op +5(t',t ))|p^ 

= W r (k)e 4k Bp+ ^ kABk+i ^(p' - p - hAk) (Dll) 



we omitted the (t',to) to alleviate the notations. The 
transition amplitude (D7l becomes: 



(D9) 



b (1 Hr,t) = in I fn _^* 2 I dt^fe,r|C/ ( i ,i )|6,p + ^',i )k)^(k) e ik -( s (*'^)P+^ t '*°)- r ») 



(27Tft)3/2(27r) 3 / 2 

xe i 2^ kA ( t '' t °) B ( t ''* ) k e^ 6a ^' _to ^F(t')e~ i 3( t '~ t ^) 2 ~ iw 2i( t '- t o)+ i vo^ a; p|\]/f£ )\ (D12) 
The evolution of the ABCDt; parameters in a constant gravitational field is simple: 

A(t',t ) = l B(t',t ) = t -^ at',t ) = l:g(t'-to) 2 (D13) 

m 2 

One can extend the computation to include the effect of interactions by taking into account an effective lensing 
effect in the ABCD matrices. This approach will be developed elsewhere. We introduc e the Fourier transform 
F(oj) = J JLF(t)e- luJt of the slowly varying envelope F(t'). The phase in the integral (ID 12b is a second-order 

polynomial in t': 

tp(t') = i(k • g - r)t' 2 + [uo ba + ^ + ^- - wax + r U - k • g t + u]t' + <p' (D14) 
2 to 2to 

where ip' is a constant phase term. In order to maximize the transition amplitude, the chirp rate r should be adjusted 
to cancel the quadratic variation of the phase, which yields as anticipated in Section [T] 

r = kg (D15) 

The time t r , at which the frequency ramp begins, selects the velocity class of the atoms which undergo the transition. 
Indeed, the momentum p of these atoms satisfies: 

k • p hk 2 

\uj ba + — £ + - u 21 +k-g(t r -t )\ <Auj (D16) 

TO 2to 
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where Aw is the spectral width of the time envelope F(t'). As a consistency check, we see that this condition 
reproduces the resonance condition ([2| for t r = to. The next step in the computation of the amplitude (D12l is to 
consider that the matter wave is out of the interaction zone at the initial and final times. This is legitimate, since we 
are in fact interested in computing a scattering amplitude. This simplification allows us to extend the bounds of the 
time integral to infinity, which yields a Dirac distribution: 



d 3 pd 3 k 



xiy(k)e- ,k - r »(a,p|*(t )) / duF(w)8[w + w ba + — - + 



(27rft) 3 / 2 (27r) 3 / 2 
k p Kk 2 



in 2m 



(b,r\U (t,t )\b,p + hk) 



^21 + k • g (t r - t ) 



(D17) 



The Bragg resonance condition selects the Fourier component of adequate frequency in the temporal envelope. To 
alleviate the notations, we note ws(p,k) the frequency selected by the Bragg condition: 



k • p Kk 2 

we(p, k) = uj ba H h w 2 i + k • g (t r - t ) 

m 2m 



(D18) 



To simplify the computation, it is useful to assume that the spectrum F(u>) is broad enough to override dispersion 
effects of the laser wave. In other words: 



ws(P)k) ~ w B (p z ,k ) 

This is legitimate if the spectral width Aw of the pulse F(t) verifies: 



(D19) 



m m 



(D20) 



Within these conditions, the Dirac integral leaves the amplitude: 

d 3 pd 3 k 



b^(r,t)=mV2^e i ^ J (^t/^wa ( h > r l^(t, t )\b, p + hk)W(k)e- ik - r ™F(u B (p z , ko))(a,p^(f )) (D21) 

In order to see how the curvature of the light beam is imprinted onto the atomic beam, we need to perform the 
integration over W(k). To compute this spatial Fourier transform, we go back to the expression of the potential (p5|: 



Vja(r,t) = ~hnU+ 2 (r-r w )e 



J2k (z-z m ) 



F{t)e 



-i(ui 2 l+r(t-t r ))t 



C.C. 



(D22) 



where we have used the relation between the Gaussian modes U + (r) = U *(r). The spatial function inside the 
potential is defined by: 



U + (v) 



1 - 2iz/b 



exp 



x 2 + y 2 



1 — 2iz/b Wq 



(D23) 



with the confocal parameter b — how 2 ,. It will be useful to introduce its transverse Fourier transform: 



U o( r ) = ^ / dk x dk y exp 



(k 2 x + k 2 )w 2 t 
■ — 4 y ° (l-2iz/b) 



(D24) 



It is convenient to introduce the Lorentzian function 



L+(z) = 



1 - 2iz/b 

The transverse Fourier transform of VF(r) can then be expressed as: 



(D25) 



W(r) = L+ (z)exp 



d 2 k± ( Wq 

W) \^ 



L + (z) exp 



{kl + k 2 ) w 2 



(1 - Hz jb) 



e i2k z \ e ik±-r± 



(D26) 
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From this last expression we infer: 
The propagation of a plane wave in a gravitational field yields 



4 2 1 7 7 



(D27) 



(r|[/ (Mo)|p) 



.iP' r P i ' s 'ci(*ito) 



(D28) 



with Sci(t,to) classical action between a trajectory of initial momentum p, final position r and duration t — t . The 
corresponding expression can be recast as: 



(r|[/ (Mo)|p) 



1 



'2irh 



_ e -img'{t-t ) /6/lg- 2^p-g(t-*o) e i(p-mg(t-t ))-r e -ip (i-t )/2ro 



(D29) 



This gives an expression for the matrix element: 



(b,r\U (t,t )\b,p + hk) = 



1 



V2tt/i 



D -img 2 (t-t f/Hh -jL p.g(t-t ) 2 _ * (p_ mg (t_t )). r ik-(r-ig(t-io) 2 -^(t-to)-r w ) 



xe 



-iu> b (t-t ) -ifr,^ 2 /2m(t-t ) -ip 2 (t-t )/2m -ihk ± 2 (t-t ) /2m 



(D30) 



The momentum distribution W(k) peaked around the value k = 2fcou z has a width Ak <C fco, and one can verify on 
(D26l that its longitudinal width Ak z is much narrower than the transverse ones Ak x ,Ak y . The correction to the 
longitudinal recoil Hk 2 /2m when k va ries i n the width of IF(k) is thus typically much smaller than the transverse 
recoil. The term Hk 2 /2m in equation (D30) will therefore be approximated by 2hko/m. By summing up the Fourier 
modes, we will recover for the atomic wave the transverse Fourier profile of the Gaussian laser wave ( D26 1 up to a 
translation. We note r c (t) = r w + ^(t — to) + |g(t — to) 2 the point associated with a classical motion in the gravity 
field from r w . The role played by the position associated with the classical movement and the action phase pre-factor 
are indeed a consequence of the ABCD^ theorem [TBI . Gathering all the terms of (D21 ) dependent on the wave- vector 
k, and using relation (D27lwe can perform the integration on the wavevector k along: 



d 2 k ± 

(27T) 



dk 



(27T) 1 / 



\j- 2 W{k)e 



ik„{z—zc) 



e ik_L-(rx-r C j.) e -i^ 



L+{z - 2c (*)) e i2 M*-*c) e -^t*-to)^ 



L + (z - Zj (t))c i2k "^ z ~ zc ~ 2hko(t ~ to)/m] 

l+(z - z c (t)y 2k ^- zo ^-^ t - t ^u+ 2 ( 



d 2 k ± 



exp 



exp 



4 2 



(1 - 2i{z - z c )/b) 



-(t-io)gikx-(ri-r C j.(t)) 



(i 



(27T) 

. . 2fi,ko , , 
r-r c (t) "(t-to) 



2i / 
k wi 



2hk (t - t ) 



!ki-(ri-r c 



(I 



The momentum acquired during the Raman process is reflected in the factor e l2fc o[z-^c(t) ^(t-to)]^ xhe translation 
r — Tc(t) — — 1 ) accounts for the classical motion in the gravitational field and the momentum acquired during 

the Raman process. Inserting this result in the equation (D21l: 

b^(r,t) — iQy/2ne i2ko ^ z ~ zco ^~ 2hk ° ( ' t ~ to ^ n ^e iv ' 



d :i i 



{2-Khfl 2 
x L+{z-z c {t))U+ 2 



e -irng 2 (t-t ) 3 /m e - ^p-g(t-t ) 2 e -ip 2 (t-t ) /2m e jr(p-mg(t-t ))-r e -iuj b t 



r - r C o 



( t ) _ 2_*o (* _ ^ F( WB (p z , ko ))( a , p|$( io )) 



(D32) 



If the atomic wave-packet is sufficiently narrow, the Gaussian modes L + (z — zc{t))U2 (r — rc(i) — ^p^ - ^o)), 
which depend on the momentum p through rc(t), are approximately constant on the width of the distribution 
F(ojb(Pz, ko))(a, Pl^^o)) centered on po- We can then pull those functions out of the momentum integral. The 
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phase factor in the momentum integral correspond to the propagation of plane waves in a gravitational field. One 
can thus interpret the momentum integral as the propagation of the filtered wave-packet in the gravitational field: 



b w (r,t) = iQ^^ k o^-^co(t)-2hk (t-t )/m] e i V ' oL +^ z _ Zco (tyu+2 ^ 



r-rco(t) -(t-t ) 



x(b,r\U (t,t ) 



d :i i 



(2irh) 3 / 2 



F(w B (p„ ko))(o,p|*(to)) 



\b,p) 



with rco(t) — r w + — (t — to) + §g(t — to) 2 - To first order in the field, the curvature of the Gaussian Raman wave is 
transferred in a controlled way to the atomic wave through the terms L + (z — zco(t))U^~ 2 (r — i"co(t) — "*K*o(t — t )). 
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